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(N : ABSTRACT 
> ■ 

^ • Images obtained by drift-scanning with a stationary telescope are affected 

^ ■ by the declination-dependent curvature of star trails. The image displacement 

• to curvature and drift rate variation increases with the angular field of view 
^ ' and can lead to significant loss of resolution with modern large-format CCD 
^ ■ arrays. We show that these effects can be essentially eliminated by means of 
1^ • an optical corrector design in which individual lenses are tilted and decentered. 

I ' A specific example is presented, of a four-element corrector designed for the 

2 ' Large-Zenith Telescope. The design reduces curvature errors to less than 0.074" 

^ ■ over a 10' x 20' field of view centered at 49° declination. By changing the 

• positions and tilts of the lenses, the same design can also be used for any field 

• centers between 0° and ±49°. 

.Si 

Subject headings: telescopes, instrumentation: miscellaneous, surveys 



1. Introduction 

In recent years the technique of time-delay integration (TDI, also known as drift- 
scanning) has become well established in astronomy ( |McGraw, Angel, fc Sargent 1980 



Hall fc Mackay 1984|) . For this technique, the telescope remains stationary, or tracks at a 



non-sidereal rate, resulting in a steady drift on images across a CCD detector (which is 
aligned so that the drift direction coincides with the orientation of the CCD columns). The 
CCD is scanned continuously at a rate which moves the photon-produced charge along the 
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columns at the same speed as the optical image. This prevents charge spread and results in 
sharp images with an integration time equal to the drift-time of the image across the CCD. 

TDl offers several benefits, compared to conventional imaging. A chief advantage 
is that variations in pixel-to-pixel sensitivity are greatly reduced because response is 
averaged over all CCD pixels in each column. The required fiat field correction is essentially 
one-dimensional. This results in as much as a factor of 45 reduction in background 
variations, after fiat field correction, for a 2048 x 2048 pixel CCD. Because background 
variations are often a limiting factor in photometry of faint or low-surface-brightness 
objects, even conventional telescopes are sometimes operated in drift-scan mode in order to 
take advantage of this and other benefits (, [Schneider et al. 19941 ). 

The TDl technique is particularly important for telescopes which have no tracking 
capability. Liquid-mirror telescopes (LMTs) employ rotating primary mirrors surfaced with 
a liquid metal such as mercury. Because the rotation axis must be vertical, these telescopes 
are best suited for observations near the zenith, although it is possible in principle to 
observe at large zenith angles using specially-designed correctors ( [Richardson and Morbey 



1987[ , [Borra, Moretto fc Wang 1995[ ). By operating their CCD cameras in TDl mode. 



such zenith-pointing telescopes can survey large areas of sky as the Earth rotates. For 
example, the UBC-NASA Multiband Survey (UNMS, [Hickson fc Mulrooney 19981) employs 
a 2048 X 2048 pixel CCD camera at the prime focus of the 3-meter LMT of the NASA 
Orbital Debris Observatory ([Potter fc Mulrooney 19971) . This survey covers 20 square 



degrees and reaches ttiab — 20.4 in a single TDl scan through medium-band filters. 

For TDl observations, it is essential that the images of all objects drift at the same 
rate, on parallel linear tracks, or image spread will result. For this reason it is important 
that telescopes used for TDl observations be equipped with correctors that reduce the 
pincushion or barrel distortion that is usually present. For a prime- focus configuration, this 
normally requires a corrector having at least four optical elements. 

After removal of telescope distortion there remains, however, a more fundamental 
difficulty. With the exception of objects on the celestial equator, the image tracks are not 
straight, but concave, in the direction of the nearest celestial pole. Furthermore, the linear 
rate of image motion is a function of declination, so there is no universal CCD scan rate 
which will prevent image spread. These effects can lead to significant image degradation for 
telescopes at moderate latitudes employing large-format CCD cameras. 

An elegant solution to this problem has been realized with the Great Circle Camera 
( [Zaritsky, Shectman fc Bredthauer 1996| ) which uses a combination of rotation and 
motion in declination to trace a great circle on the sky. As a result, the nonlinear effects 
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are minimized. While the Great Circle Camera does reduce curvature effects, it has a 
fundamental limitation. Because of the physical motion required of the camera, the length 
of the great circle track is limited. This restriction is particularly severe for zenith-pointing 
telescopes which cannot move to follow a great-circle track. TDI observations would need 
to be interrupted frequently while the camera position is reset to follow different great-circle 
arcs. 

In this paper we propose a new solution to the problems of star-trail curvature. By 
offsetting and tilting the lenses of a suitably-designed corrector, one can introduce an 
asymmetric field distortion which compensates for the nonlinear motion of the images. As 
a result, the images track in parallel straight lines at a common constant rate with high 
accuracy. This effect can be achieved without introducing any additional optics in the 
light path. Moreover, the same corrector can serve telescopes at different latitudes just by 
repositioning the optical elements. 

The paper is organized as follows: In Sec. 2 we analyze the geometry of star-trail 
curvature and derive formulae for the resulting image displacements. In Sec. 3 we describe 
an example of a prime-focus corrector which compensates for these effects. The performance 
of the corrector is examined in Sec. 4 followed by our conclusions. 



2. Analysis of Star- Trail Curvature 

Star-trail curvature effects have been discussed previously by ( Gibson fc Hickson 1992 



Stone et al. 199(^ and [Cabanac 1998| ) who give approximate formulate for the distortions. 



Our aim here is to present a simple, exact, derivation of the relevant relations as well as 
useful approximate formulae. The geometrical picture is illustrated in Figure 1. The circle 
represents the celestial sphere, which for convenience we take to have unit radius. Its center 
is the origin of a three-dimensional Cartesian coordinate system {x,y,z). The z-axis is 
aligned with the North Celestial Pole (NCP) hence the x-y plane intersects the sphere at 
the celestial equator. For clarity, the y axis, which is perpendicular to the page, is not 
shown in the figure. Now consider a telescope whose axis lies in the x-z plane at an angle 
6o to the celestial equator. Sq is therefore the declination of the center of the telescope field 
of view. For a zenith-pointing telescope, 5o is equal to the latitude of the observatory. In 
the telescope's focal plane, perpendicular to the optical axis, we set up a two-dimensional 
Cartesian coordinate system {X, Y) , oriented so that the F-axis is parallel to the y-axis. 
The optical axis passes through the origin of the X-Y coordinate system, so that X and 
Y give the North and East position of the image with respect to the field center. It is 
convenient to work with the projection of the focal plane, through the telescope optics. 
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back onto the celestial sphere. Figure 1 shows the projected focal plane which is tangent to 
the celestial sphere at the point where it is intersected by the optical axis. Since the sphere 
has unit radius, the values of X and Y can be converted to physical distances in the focal 
plane by multiplying by the telescope's effective focal length F. 

Now consider the sidereal motion of a star at arbitrary declination 6. As the Earth 
rotates, the line of sight to the star sweeps out a cone, indicated in Figure 1 by the '*' 
symbol. The intersection of this cone with the tangent plane defines the track of the star's 
image on the focal plane. This track is by definition a conic section, which is an ellipse, 
parabola or hyperbola depending on whether the declination 6 is greater than, equal to, or 
less than the co-angle ir — Sq (for a zenith telescope, the co-latitude). 



The equation of this track is easily derived. The equation of the cone is 

z'^ cot^ 6 , 



2 I 2 

X + y 



while the plane is defined by 



X = —X sin Sq + cos Sq 
z = X cos 6o + sin 6o ■ 



Eliminating x and z, and setting y = Y, we obtain 

(sin^ 6 - cos^ 6o)X'^ - 2 sin 6o cos 6oX + (F^ + 1) sin^ 6 - sin^ 6o = 



(2) 



(3) 



The northerly displacement X of the image as a function of East- West position Y is given 
by the solution 



X 



sin 6o cos 6o 
sin^ 6 — cos^ 5n 



1 



1 - 



(sin^ 6 - cos^ 6o){[l + Y^] sin^ 6 - sin^ 



sin 6o cos^ 5o 



1/2 > 



This equation has a simple Taylor-series expansion: 

y2 



X = tan(5 - 6o) 



2 cot (5 



+ 



(4) 



(5) 



From this we see that the local radius of curvature R measured at the center of the 
track is 



R = 



d'X 



dY^ 



-1 

Y=Q 



cot 6 



(6) 



which in physical units is F cot S. Note that this curvature depends only on the declination 
of the object, and is independent of the pointing direction or location of the telescope. 
Equation (^) may be written 

X = Xo + AX , (7) 
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where Xq = tan{S — Sq) denotes the X-coordinate of the center of the track {Y = 0), 
and AX measures the departure of the track from hnearity, in the North-South direction. 
Except at very large dechnations, the radius of curvature of the trails is nearly constant 
across the field and we may use the value at the field center: 

AX ~ ^r^tan^o . (8) 

Consider now the rate at which the images move. Let be the azimuth angle of the 
object in the (x, y, z) coordinate system, referenced to the x-axis. As the Earth rotates, 
increases at a constant rate. From the definition of 0, 

Y = y = xtan0 . (9) 

Substituting this in Eqn (|^) and using Eqn (|^) to eliminate we obtain 

Y = . (10) 

sin (5o tan b + cos (5o cos 

From this we see that the rate of image motion is a function both of the declination 
of the object and the declination of the field center. For a given azimuth angle 0, objects 
which cross the field North of the field center (Xq > 0) will have smaller Y values that those 
passing south of the field center (in the Northern Hemisphere - the opposite is true in the 
Southern Hemisphere). 

Let 1^ = 1^0 + where Yq denotes the ^-coordinate of a star, having the same 
azimuth angle, whose track passes through the field center (Xq = 0). From Equations (|^), 
(0) and (|10D we obtain 

AF = -Xpyoj l + t^"^^;^^o I . (11) 
\ tan + cot oq cos J 

The factor in parenthesis typically varies very slowly across the field and may be replaced 
by its value at the field center. This gives 

Ay ~ -Xoyotan(5o . (12) 



3. Correcting Curvature Effects 

In the previous section, exact and approximate formulae were developed to describe the 
distortions produced by star-trail curvature. For most practical situations, the approximate 
formula of Eqns (§) and |T^) give more than enough accuracy. The maximum distortion can 
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be found by setting Xq and Yq equal to half the North-South and East- West extent of the 
detector, respectively. The two distortions are closely related. For a square detector, the 
East- West image smear, due to rate variations, is four times as large as the North-South 
smear due to curvature. (Because the distortion changes sign as the star crosses the field, 
the East- West image smear is twice AY.) 

As an example, consider observations of a star at 5 = 45° using a CCD camera 
with a square field of view of 20' (Xq = Yq = 0.002909 radians). Equations (g) and (0) 
give maximum distortions of 0.000004231 radians and 0.000008462 radians (0.873" and 
1.745") respectively for AX and AY. These are significant distortions which would cause 
unacceptable image degradation for TDI observations 

Fortunately, it is possible to largely eliminate these effects for zenith telescopes, located 
even at moderately-high latitudes, by means of an asymmetric corrector. The idea is to use 
a combination of decenters and tilts of the optical elements in order to introduce distortion 
that closely matches that described by Eqns (|D and (|12[), but with opposite sign. Such a 
design has been recently developed for the Large Zenith Telescope (LZT), which employs 
a 6 meter f/1.5 liquid-mercury primary mirror ( piickson et al. 1998| . The corrector was 
designed by E.H.R. to specifications provided by P.H. The required image quality was a 
50% encircled energy diameter (EED) of 0.4" or less over a 10.5' x 21' field of view (to 
match a 2048 x 4096 pixel CCD). The procedure was to first aim for a corrector with 
zero distortion. Global optimizations were employed to keep the maximum clear aperture 
relatively small, in order to reduce cost, while meeting image-quality specifications over 
the required field of view. The corrector was then reoptimized to introduce the required 
sidereal distortion. In this process, curvatures, tilts, decenters and spacings were all allowed 
to change. At first, a wedge was allowed for one element, but it was found that this could 
be eliminated, allowing the same lenses to be used at different latitudes. After optimization 
for the +49°17' latitude of the LZT, the lens curvatures and thicknessess were frozen. In 
order to illustrate the flexibility of the corrector, the design was then reoptimized for the 
+32°59' latitude of the NODO LMT, and the equator (0°), allowing only the tilts, decenters 
and locations of the lenses to change. 

The optical configurations, for latitudes +49°17' and -|-32°59', are illustrated in Figures 
2 and 3. The corrector employs four lenses all fabricated from the same optical material. 
The details of the corrector design are provided in Table 1, which gives parameters for the 
two asymmetric configurations and for the symmetric case. Each line of the table describes 
an optical surface, in the order in which light reaches it. Light enters parallel to the z axis of 
a Cartesian coordinate system aligned with the Surface 1 (the parabolic reflector). Columns 
(7), (8), (10) and (11) describe translations and rotations of the coordinate system in which 
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subsequent surfaces are defined. Surfaces following such a transformation are aligned on 
the local mechanical axis (z-axis) of the new coordinate system. The new mechanical axis 
remains in use until changed by another transformation. In these transformations, the 
translation is applied before the rotation. 

The column headings for the table are as follows: (1) surface number, (2) material 
that the light enters when crossing the surface, (3) radius of curvature of the surface, (4) 
the aspheric constant A (except for surface 1 which is a pure conic section), defined by the 
equation 

Az = ^ + Ar^ , (13) 

l + {l-(l + X)yV7?2}V2 

where r is the distance from the z-axis to a point on the surface and Az is the 2;- displacement 
of the surface with respect to a plane, (5) distance along the z-axis from this surface to the 
next for the 0° configuration, (6) change in the z-axis distance required for the -|-32°59' 
configuration, (7) displacement of the origin of the coordinate system at the surface, (8) 
rotation angle of the coordinate system at the surface. Columns (9), (10) and (11) are the 
same as columns (6), (7) and (8) but for the -|-49°17' configuration. 



4. Performance of the Corrector 

The corrector provides a 24 arcmin-diameter unvignetted field of view. The effective 
focal length is 10.000 m which gives an image scale of 20.63 arcsec per mm. Table 2, gives 
the 50%, 80% and 100% EEDs for various field angles and wavelengths. The median values 
of the 50% EED at 400 nm, 500 nm and 600 - 800 nm are 0.353, 0.162 and 0.203 arcsec 
respectively. These values are well below typical ground-based seeing disk diameters at 
these wavelengths, which are typically 0.6 arcsec or more at the best astronomical sites, so 
the corrector will not appreciably degrade the image. There is a small, but non-negligible, 
focus shift between the three wavelength regions. Since the the corrector is intended to be 
used with common broadband (or narrower) filters, this is not an issue for image quality. 

Table 3 summarizes the distortion characteristics of the corrector. For the field angles 
specificed in Columns (1) and (2), Columns (3) and (4) give the target and achieved image 
displacements. The difference in position, which corresponds to the distortion error, is 
listed in Column (5). It can be seen from the table that the distortion errors are typically a 
few tens of milliarcsec, the maximum error being 0.074". Since these residuals are much less 
than the seeing diameter, the corrector effectively eliminates star-trail curvature and rate 
differentials as a source of image degradation. 
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Considering the 6 meter diameter and fast focal ratio of the primary mirror, the 
corrector is very compact having an overall length, from first element to focal plane, of 50.5 
cm. The largest lens, which has spherical surfaces, has a diameter of only 34 cm. This is 
quite small for a 6-meter f/1.5 telescope, which helps reduce the cost of the corrector. The 
three smaller lenses have diameters less than 15 cm. The rear side of each is a low-order 
aspheric surface. 

5. Summciry and Discussion 

We have presented an analysis of star-trail curvature effects. If uncorrected, these 
effects can result in substantial image degradation in drift-scan images obtained with 
large-format CCDs. A new technique has been described in which the curvature effects arc 
compensated by means of a corrector lens employing decentered and tilted elements. As an 
example, we have discussed a compensated prime-focus corrector designed for the 6- meter 
LZT. 

The corrector requires no additional optical components other than the four elements 
normally needed to remove telescope aberrations, including distortion. By varying the 
positions of the the elements, the corrector can be used with a zenith-pointing telescope at 
any latitude up to at least ±50°. It provides a 24' unvignetted field with median 50% EED 
of 0.2" and a maximum distortion error of order 0.07". 

An interesting question is how much image quality is sacrificed in order to achieve 
the distortion correction. The instantaneous image quality of the 49 deg configuration is 
slightly worse than that of a similar corrector optimized for zero distortion. However the 
difference is much less than the improvement in integrated image quality provided by the 
distortion correction. For example, reoptimizing the corrector for zero-distortion, allowing 
the lens shapes to change as well as the separations, results in only a 25% improvement in 
the worst-case RMS image spot diameters. 

Can this design be extended to larger fields of view, and how would the result compare 
to that of the Great-Circle camera? While we have not explicitly investigated this question, 
during the course of our work an earlier design was made for a corrector which has a 30' 
diameter field of view (to accommodate a 4096 x 4096 pixel CCD) and assumes a 5-m f/1.8 
primary mirror. The image quality at 49° is comparable to that of the present design. 
Because of the smaller primary mirror, a strict comparison with the present design is not 
possible, but it does show that wider-field designs are possible. 

In this context, it should be pointed out that even the Great-Circle camera is not entirely 
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distortion-free. The image displacements can be obtained from Equations (y) and (11) by 
setting do = and taking 6 to be the field angle X. This gives AX = tan 5/2 ^ XY'^/2 
and Ay = 0. So while the great-circle camera is free from rate-variation, there is a small 
amount of star-trail curvature. This results in an image smear which increases in proportion 
to the cube of the field diameter and becomes substantial for field sizes of order one degree 
(the image displacement is 0.548" for a 1° x 1° field). 

While the corrector design presented here was specifically developed for a zenith- 
pointing telescope, the same optical design could also be used with a conventional telescope. 
Because the curvature effects depend only on declination, compensation could be introduced 
for any field center by means of mechanical actuators that would adjust the decenters and 
tilts of the individual elements in the corrector. The corrector would be reconfigured for 
each pointing of the telescope, and then held fixed during the integration. This would be 
the case even if the scanning is done at a non-sidereal rate, so long as the declination of the 
field center remains constant during the exposure. 



We thank the referee. Dr. Stephen Shectman, for helpful comments. Research support 
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of Canada (NSERC). The LZT is funded by NSERC Collaborative Project Grant no. 
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Fig. 1. — Geometry of star-trail curvature. 

Fig. 2. — The LZT corrector configured for observations at +49°17' declination 
Fig. 3. — The LZT corrector configured for observations at -|-32°59' declination 
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Table 1. Corrector Specifications 
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Table 2. Encircled Energy Diameters for +49° 17' 
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Table 3. Image Displacements for +49°17' 
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0.565 


0.554 


0.023 


55 


11 


600 


0.283 


0.283 


0.565 


0.526 


0.039 


11 


11 


800 


0.283 


0.240 


0.565 


0.513 


0.066 


0.044 


0.088 


400 


0.283 


0.259 


-0.565 


-0.605 


0.047 


55 


55 


600 


0.283 


0.259 


-0.565 


-0.569 


0.024 


55 


55 


800 


0.283 


0.259 


-0.565 


-0.553 


0.027 
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